Theory of inter-edge superexchange in zigzag edge magnetism 
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A graphene nanoribbon with zigzag edges has a gapped magnetic ground state with an antiferromagnetic 
inter-edge superexchange interaction. We present a theory based on asymptotic properties of the Dirac-model 
ribbon wavefunction which predicts W^^ and W^' ribbon-width dependencies for the superexchange interac- 
tion strength and the charge gap respectively. We find that, unlike the case of conventional atomic scale superex- 
change, opposite spin-orientations on opposite edges of the ribbon are favored by both kinetic and interaction 
energies. 



Introduction — Motivated by the seminal theoretical work of 
Kobayashi, Fujita, Wakabayashi and collaborators|[ll|2,[3ll3], 
and by progress in gra phene p reparation||5i1, researchers have 
recently S 11 S |1 E [iil III H H El 111 reexamined 
the intriguing physics of edge magnetism in zig-zag termi- 
nated graphene nanoribbons from a number of different points 
of view. The magnetic state is a consequence of the nearly- 
flat subbands which occur at the Fermi level in a neutral 
zig-zag ribbon, and of the orbital character of the wavefunc- 
tions associated with these bands. Although there is still no 
conclusive 1 18] experimental evidence that a one-dimensional 
(ID) magnetic state occurs in ideal zigzag ribbons, the the- 
oretically predicted state seems likely given that quite differ- 
ent electronic structure theories (from crude Hubbard models 
to elaborate ab initio DFT calculations) yield consistent! 14] 
predictions and that there are at present no other ideas on 
how the unusual flat bands could be accommodated in the 
many-electron state. Present ribbons are far from ideal, how- 
ever, and the main obstacle to realizing this paradigmatic ex- 
ample of d'^ magnetismi 19, 20] may lie in furthering recent 
progress fl2ll l22i ,23 1 toward chemistry and defect control at 
the edge. 

Mean-field-theory calculations predict that the ground state 
of a zigzag ribbon has unusuafly stiff parallel spin-alignment 
along each edgefis'l and antiferromagnetic| Ts'l inter-edge su- 
perexchange interactions. In this Letter we present a mostly 
analytic theory of the ribbon width W dependence of the im- 
portant inter-edge superexchange interaction. Our theory re- 
lies on the properties of large W solutions of the continuum 
model approximation 101 for zigzag edges. We predict that in- 
teredge interactions can have a substantial influence on the 
properties of these unusual ID magnets. 

Ribbon edge-state bands — The TT-orbital tight-binding model 
for a finite width ribbon yields a number of one-dimensional 
Bloch bands proportional to the ribbon width lfl6ll . To a high 
degree of accuracy, zigzag magnetism is a rearrangement of 
only lfl4ll the highest occupied (|^ — )) and lowest unoccupied 
{\k +)) ribbon bands, whose transverse wavefunctions are re- 
spectively odd and even functions of carbon atom sites across 
the ribbon. Since the exchange physics which favors mag- 
netism is local, it is revealed most clearly by forming states in 
this Hilbert space which are localized as far as possible at one 



edge or the other: 

\kL) = ^i\k-) 



-\k+)), \kR) = ^{\k+)-\k-))(l) 



where we have chosen the band transverse wavefunction am- 
plitudes to be positive at the left most atom. These L(eft), 
R(ight) basis states can be expanded in terms of amplitudes 
on atoms in the ribbon unit cell: \kL) — Liii\kl)g, \kR) — 
Rki \kl)g with sums over sites / within the unit cell implied. 
It is readily verified that L^i and /?^/// are strictly locatized 
on opposite sublattices so that Li,iRyii = 12411 . In this 
representation, the tt -band tight-binding model Hamiltonian 
Htb {k) = t (k) where is a Pauli matrix and the left-right 
tunneling amplitude t {k) ~ —Ejb (^) ^ ^tb (^)- The bands 
in the Brillouin zone —n/a<k< n/a have periodicity lit /a 
and have inversion symmetry so we can restrict our attention 
to the interval < k < n/a. Zigzag edge magnetism fol- 
lows from the following tight-binding model property 10, [H]. 
For 2k/3 + qea < \k\, where qe — i /W — 2/ \/3aN, the states 
\kL) and \kR) are exponentially localized near their respec- 
tive edges and the left-right hopping amplitude t [k) decreases 
rapidly with increasing W . (Here a = 2A6A is the lattice con- 
stant and is the number of atom pairs in the ribbon unit cell.) 
Over this region of wavevector \kL) and \kR) are proper edge 
states. 

Hubbard model mean-field theory — In order to explain our 
theory of the superexchange interaction we briefly summarize 
the mean-field theory of the magnetic state, which is partic- 
ularly simple in the Hubbard model case. It is instructive 
to contrast two different collinear magnetic solutions of the 
mean-field equations, an antiferromagnetic (AF) one in which 
spins have opposite orientations on opposite edges and a ferro- 
magnetic (F) one in which spins have the same orientation on 
opposite edges. The two LR basis spin-dependent mean-field 
Hamiltonians are given (up to a common constant) by: 
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where a = +/— for t / i spin, the A:-dependence of t and the 
self-consistent exchange potentials is implicit, and 



{k) = uY,Lh{mi)^P, 
I 



A^{k) = U'£Lli{m,)p. (3) 
/ 
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In Eq.([3]l mi = {n^ — ni\^ /2 is the site-dependent spin- 
density and n/fj is the spin-dependent mean occupation num- 
ber at site /. (Both solutions have «/ = («/|+n/j) = 1, a 
convenient property of neutral ribbons which can be traced 
to particle-hole symmetry in the paramagnetic bands.) These 
self-consistent solutions are illustrated for a. N — 20 zigzag 
ribbon in Fig.( [T]). Both the F (even m/) and AF (odd m/) 
state solutions are self-consistent. The mean-field equations 
are closed by evaluating ot/ using 



±(2;r/3 +q) and small q the expression 

f(±(^ + 2;r/3)) = 
where z satisfies 



zWcoth(zW) 
zW/tan(Wz) =^W, 



q> qe 

q < qe- 



(7) 
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In the continuum model 
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("^')af = 2^ dk{Ll,^Rl,) P{k) (4) /^f^ ^ ^ 2z (cosh(2z3;) - 1) / (sinh (2Wz) - 2Wz) (9) 

('«/>F = 



^ I dk (rI 



where 



\AF 



P{k). 



(5) 



(6) 



is the degree of left-right edge polarization of the AF mean- 
field states. In the AF case, local spin-polarization follows 
from the opposite left-right polarizations of f and | states 
whereas in the F case the left-right polarization vanishes and 
spin-polarization follows from double occupation of f bands 
for kc < \k\ < n/a where kc is the wavevector at which the 
|,+ band and the J,, — band cross as illustrated in Fig. ([T]). 
A'*^ and A^ are nearly (but not quite!) identical because their 
m/'s differ mostly in the middle of the ribbon which has httle 
influence on the edge states. 
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FIG. 1: (Color online) Hubbard model mean-field calculations 
for 7() = 2.6eV and U = 2eV for a zigzag ribbon with N = 20. 
Left panel: Mean-field energy bands for the AF state: E'^^ (k) = 
1 /2 

± ^A'*^ (fe)2 -|- (fe)^ for both spins with the low-energy states of 

one spin concentrated on one side of the ribbon and the low-energy 
states of the opposite spin concentrated on the opposite side. t{k), 
the left-right hopping parameter, is quite insignificant for this ribbon 
width in the edge state region. A"^^ (k) is dominant in the edge state 
region because of large local spin-polarizations. Right panel: Mean- 
field energy bands for the F state: (k) = ctA'" (k) ± \t {k)\. Note 
that A'*^ and A'^ are nearly identical. The bands are periodic with 
periodicity 2n/a and inversion symmetric. 



Ribbon width scalins rules — From solutions of the graphene 
continuum model 10, [l6ll we obtain for the region near k — 



for^^. <k <n/a and 

Rl{y) = 4zsin {zy) / {2Wz - sin {2Wz) ) . 



(10) 



for < A: < ^c- The left centered functions can be obtained 
through the symmetry relation Lj^i — RkiN-i- It follows that 



and that 



t{k)^^t{qW) 



(11) 



(12) 



where the functions X^(^) t{x) are impUcitly defined by 
the above equations. We have verified that these scaling rela- 
tions apply accurately even in quite narrow ribbons. 
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FIG. 2: (Color online) Left panel: Dependence of A'^'^ and /q/ on the 
scaled coordinate q = qW and the corresponding N = 60 AF state 
quasiparticle bands. Note that the self-consistently calculated A'^'^ 
approaches a well defined function at large N. The positions kc and 
k* are respectively the values of k at which A'^'' = t and the band gap 
minimum occurs. Right Panel: Scaling collapse^of antiferromag- 
netic state self-consistent left-right polarization P and symmetric- 
antisymmetric polarization T represented in the scaled coordinate q. 
Note that P^ + = Ihy definition. 



From Eq.lO we see that the Hubbard model exchange po- 
tentials depend on local spin polarizations (m/) which are 
large only close to the edge and approach a well defined 
limit already for quite narrow ribbons; the form of the spin- 
polarization near each edge is a single-edge property unrelated 
to interedge interactions. From this observation and the above 
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scaling relations for the zigzag edge states, we propose the 
following scaling rule for the form of the exchange potential 



Since both A (k) and t (k) depend only on qW it follows from 
Eq.(|6]l that P{k) also depends only on qW. We have verified 
numerically that this scaling relationship holds accurately for 
sufficiently wide ribbons as illustrated in Fig. (El. 
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FIG. 3: (Color online) fc-resolved contributions to the kinetic (ej) 
(upper panel) and exchange (lower panel) energies (ex) of the F 
and AF states and the corresponding F — AF differences as a func- 
tion of the scaled momentum coordinate q = qW. £j and e^'^ are 

the integrands in the kinetic energy expression Eq.(|15ll, and e^^^^^ 
is the corresponding quantity for the exchange interaction energy. 
The discontinuities in the ferro case are due to the crossing between 
the majority-spin symmetric and minority-spin antisymmetric bands 
at kc, indicated by a thin vertical line. Although the kinetic ener- 
gies roughly double the interaction energies at most k values, the 
exchange contribution to superexchange is much larger because of 
weaker cancelation between \k\ < kc and \k\ > kc regions. 



Inter-edge interaction — The strength of the superexchange 
interaction which determines the alignment between magneti- 
zation directions on opposite edges is given by the total energy 
difference between AF and F solutions. Because the electro- 
static Hartree energies of both states are identical, the energy 
difference per edge carbon atom A/s can be separated into 
band (kinetic) and exchange energy contributions: 

7F 



AE = 



A£v 



(14) 



The difference of kinetic energies between AF and F solutions 
are determined by contributions from occupied edge band 
states: 



j^AF ^ I dkt{k)T{k) 



la 
n 
2a 
n 



dk t (k) 



where T{k) = (1 — is the symmetric-antisymmetric 

polarization of AF states. The F as well as non-interacting 
band eigenstates have r(^) = 1. In the non-interacting ground 
state the lower energy state is fully occupied and the total en- 
ergy contains all the band energy. Both AF and F states sac- 
rifice band energy contributions in the region \k\ > 2n/3a in 
order to gain interaction energy. In the ferromagnetic case 
the band energy gain is sacrificed completely for \k\ > kc, the 
wavevector at which A = f . At larger values of \k\ both bond- 
ing and antibonding states are occupied for one spin and both 
are empty for the other spin. There is therefore an abrupt sep- 
aration at = kc between wavevectors which contribute to 
band energy and regions which contribute to the exchange en- 
ergy, discussed below. In the AF case, on the other hand, T 
crosses smoothly as a function of scaled wavevector q = qW 
from the kinetic energy contributing regime at small \k\ to 
the exchange energy contributing regime at large The in- 
teredge interaction is due to this difference. Using the scaling 
properties of t and A the kinetic energy contribution to the 
difference can be written as an integral over q: 



AT 



2a Yo 



dqt{q) [T{q)-\]- 



dq t{q) T{q) 



(16) 

The integrals in Eq.( [T6l) converge at —°o because T{q) ap- 
proaches 1 rapidly and at °o because T{q) approaches 
rapidly. The contribution from q < qc is negative while the 
contribution from q > qc is positive. Substantial cancelation 
leads to a small kinetic energy contribution to AE. 
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FIG. 4: (Color online) The difference between the F and AF states in 
total energy (A£), exchange energy (AEx) and kinetic energy {AT) 
plotted on a linear scale (left) and a logarithmic scale (right). The 
total energy difference follows a W^^ decay law at large W and 
is dominated by exchange energy contribution. The kinetic energy 
contribution is substantially smaller and the asymptotic decay law 
develops only for sufficiently large ribbon width. The total energy 
AE was fitted with 2.7/(W^ + 280), the exchange energy AEx with 
2.1/(W^ + 100) and the fitting for the kinetic contribution was ob- 
tained from the difference between both resulting in 0.6/W^ (repre- 
sented with x20 magnification), all terms given in eVA^^ units. We 
used a I2K A:-point sampling density in the ID Brillouin zone. 



The exchange energy integrands satisfy the scaling rela- 
tions similar to the kinetic terms and we can write the ex- 
change energy difference as 



4 



+ ^-{^^^{K/a)-^'{K/a)) (17) 

The first two terms are similar to the corresponding band energy contributions, with the discontinuity at qc again due to the 
band crossing in the ferromagnetic state. We write this contribution to the superexchange interaction as Jx/W'^. An additional 
contribution appears because A'^^ and A^ are not quite identical f or q^°°. 



A£v 



dqA^'{q)[P{q)-l] 



In the Hubbard model we can relate the asymptotic difference 
in A to the difference in spin polarization on the edge atom: 

5A = A-*^ (tt/o) - A^ {n/a) = U {{m,dge)Ap ~ {medge)p)- 

LabeUng the leftmost site as site 1 , noting that R^^ = and 
recalling the definitions of (m/) in Eqs. ( |4|5t we find that 



5A 



a 



dkLl^P{k) 



2n Jkc 



%/a 



dkLl,{\-P{k)) 



(18) 



Adding the three contributions, the total superexchange inter- 
action is 



AE{W)^W-'{JK + Jx+Jsn 



(19) 



For 7o = 2.6eV and on site repulsion U — l.QeV that results 
in band gaps similar to LDA il4ll . we find that the kinetic and 
exchange contributions to the interaction are 
andJx+J5„, = 2.leV-A^. (See Fig.(|4li.) Separately 75„, ~ 
1.1 SeV -A^ implying that the interaction contribution is com- 
posed in approximately equal measures of contributions from 
q near q^ and contributions far in the edge regime. 
Conclusions and discussions — Our analysis shows that the 
antiferromagnetic inter-edge superexchange interaction in 
magnetic zigzag nanoribbons is the sum of three contributions 
(band-energy, exchange-energy, and edge spin-polarization) 
all of which arise from a region of the ribbons's ID Brillouin- 
zone which is centered on \k\ — In /3a and scales in width as 
1 /W. Unlike the familiar case of atomic-scale superexchange 
interactions, in which antiferromagnetic spin-arrangements 
lower the kinetic energy at a cost in interaction energy, all 
three contributions have the same sign - with the kinetic con- 
tribution being substantially smaller in magnitude. Our con- 
clusions rest primarily on analytic properties of continuum 
model solutions to the ;r-band model for zig-zag nanoribbons 
and depend on the particle-hole symmetry of graphene's con- 
duction and valence bands. We have demonstrated numeri- 
cally that the continuum model predictions are accurate, even 
in narrow nanoribbons. Although some details of our analysis 
depend on the simplified Hubbard model we use, we expect 
that the scaling properties of states near \k\ = 2n/3a to be 
general and that our qualitative conclusions will apply to any 
mean-field-theory treatment of zigzag ribbon magnetism. 

Collective spin-behavior is expected lll3ll to be impor- 
tant in zigzag ribbon magnets, even though they are one- 



dimensional, because of the exceptionally strong exchange 
interactions along each edge. Assuming that the magnetic 
anisotropy (which is expected to be weak) is important only 
at low-temperatures, the correlation length ^ along an isolated 
zigzag edge is estimatedl 13] to be ^ 3000A/r[/r]. Since the 
interedge interaction arises from an interval of fc-space with 
width 1 /W, its range along the edge will be ^ W. When ^ 
is smaller than W, interedge interactions will have little influ- 
ence. For £, larger than W, the interedge interactions will help 
suppress thermal magnetization fluctuations. 
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